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Abstract 



We present a measurement of the electron spectrum from inclusive semilep- 
tonic B decay, using 5.1 fb _1 of T(45) data collected with the Belle detector. 
A high-momentum lepton tag was used to separate the semileptonic B decay 
electrons from secondary decay electrons. We obtained the branching frac- 
tion, B(B -> Xe + u) = (10.90 ± 0.12 ± 0.49)%, with minimal model depen- 
dence. From this measurement, we derive a value for the Cabibbo-Kobayashi- 
Maskawa matrix element \V cb \ = 0.0408 ± 0.0010(exp) ± 0.0025(th). 
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I. INTRODUCTION 



The inclusive semileptonic branching fraction of B decay has long been an interesting 
puzzle in heavy flavor physics. Bigi et al. |l] first pointed out that theoretical calculations 
including QCD corrections disagreed with experimental results. While most measurements 
have consistently been smaller than 11% 0-fl, theoretical expectations have been higher 
than 12%. Some theoretical analyses have been better able to accommodate a low semilep- 
tonic branching fraction by including final state mass effects in the next-to-leading order 
QCD corrections This does not necessarily solve the puzzle, however, because the 

semileptonic branching fraction is correlated with the rate for B(b — *> ccs), which depends 
on both the quark- mass ratio (m c /mb) and the renormalization scale (n/mb), where fi is 
used to renormalize the coupling constant and Wilson coefficients (c±(/i)) appear- 

ing in the non-leptonic decay rate. Current measurements of the charm multiplicity in B 
decay (n c ) |7| cannot easily accommodate the semileptonic branching fraction measured at 
the T(4S) resonance, although measurements at higher energies ||[| are somewhat more 
compatible. Further studies at the T(4S) are needed to achieve a better understanding 
of B decay and theoretical models. In addition, the semileptonic branching fraction may 
be combined with the decay lifetime to obtain the partial width, which is used to extract 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |V c b| [ID] and to probe theoretical 
models of B decay. 

For this measurement, it is essential to distinguish between primary decay leptons from 
B — > Xl + v and secondary decay leptons produced mainly through charm decay (B — > 
DX,D — > Y£~&). We have used a dilepton method introduced by ARGUS to minimize 
model dependence in this measurement. This approach requires a high momentum lepton, 
an electron or a muon, to identify a BB event and to tag the flavor of one of them. We then 
study the spectra of additional leptons in the event, extended to low momenta, selected so 
that they come mainly from decays of the other B, and separated by charge relative to the 
tagging lepton to distinguish primaries and secondaries effectively. This spectrum study is 
performed using electrons only, as their experimental identification extends to much lower 
momenta than that of muons. The numbers of electrons in our sample are determined by 
fitting the distributions of the ratio of cluster energy to track momentum {E/p) in each 
kinematic bin. The branching fraction is then obtained by normalizing to the total number 
of tag leptons, rather than the luminosity or the number of BB events. 

In this paper, we report on the branching fraction of inclusive semileptonic B decay 
(B — > Xe + u) and the CKM matrix element |V^,|. Charge conjugation is implicitly included. 
The data sample used in this analysis was collected at the T(4S) resonance with the Belle 
detector JTT]] at the KEKB asymmetric e + e~ (3.5 on 8 GeV) collider \T2\. This analysis is 
based on an integrated luminosity of 5.1 fb -1 at the T(4S) resonance, and 0.6 fb _1 continuum 
data at a center-of-mass energy of 60 MeV below the T(4S). A GEANT [O based Monte 
Carlo simulation was used to model the detector response. 



II. THE BELLE DETECTOR 

The Belle detector is configured around a 1.5 T superconducting solenoid and iron struc- 
ture surrounding the KEKB interaction region. It covers 92% of the total solid angle in the 
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T(4S) center-of-mass (CM) system. Charged-particle tracking is provided by three layers 
of double-sided silicon vertex detectors (SVD) and a 50-layer central drift chamber (CDC). 
Eighteen of the wire layers provide a small stereo angle to measure the coordinates of the 
particle trajectories in the direction of the beam (z). The transverse momentum resolution 
for charged tracks is (<7 Pt /pt) 2 = (0.0019p t ) 2 + (0.0030) 2 , where p t is in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the CDC, an aero- 
gel Cerenkov counter (ACC), and a time-of- flight scintillation counter (TOF). The dE/dx 
measurements have a resolution for lepton tracks of 6% and are useful for separating elec- 
trons from hadrons over nearly the full momentum range. The ACC and the TOF are used 
to reject charged K and proton tracks. 

The electromagnetic calorimeter (ECL) contains 8736 CsI(Tl) crystals located behind the 
hadron identification detectors, inside the solenoid. Its thickness is 16.2 radiation lengths 
over the entire tracking acceptance. The photon energy resolution is (<je/E) 2 = (0.013) 2 + 
(0.0007/E) 2 + (0.008/E 1 / 4 ) 2 , where E is in GeV. 

The h/Kl detector (KLM) located outside the coil consists of 15 sensitive layers and 14 
iron layers in the octagonal barrel region. The iron plates and associated materials provide 
a total of 4.7 nuclear interaction lengths at normal incidence. 

III. DILEPTON SELECTION 

Hadronic events were selected based on charged track information from the CDC and 
cluster information from the ECL after transformation to the CM system. We required at 
least three charged tracks, an energy sum in the calorimeter between 10% and 80% of \/s, 
a total energy sum of greater than 20% of y/s, and a total momentum balanced in the z 
direction within 2.1 GeV/c. This removed the majority of two-photon, radiative Bhabha, 
and t + t~ events where both r's decay to leptons. Remaining radiative Bhabha events were 
removed by requiring at least two large angle clusters in the ECL and that the average 
cluster energy be below 1 GeV. In order to remove higher multiplicity r + r~ events, we 
calculated the invariant mass of the particles found in each hemisphere perpendicular to 
the event thrust axis and removed the event if it fell below the r mass. Beam-gas and 
beam-wall backgrounds were removed by reconstructing the primary vertex of the event and 
requiring it to be consistent with the known location of the interaction point (IP). To suppress 
continuum, we required that the ratio R2 of second to zeroth Fox- Wolfram moments [|14j] be 
less than 0.6. 

Events containing a high-momentum "tagging lepton" to tag the B flavor and a "spec- 
trum electron" for the spectrum study were selected from the hadronic event sample. To 
select tagging electrons, we used dE/dx measurements in the CDC, the response of the 
ACC, matching between the ECL cluster and associated CDC track, the shower shape of 
the ECL cluster, and E/p, the ratio of the ECL energy to the momentum of the associated 
CDC track. Each set of measurements was translated into an electron probability. These 
probabilities were then combined into a single likelihood discriminant, L e . For the selection 
of tagging muons, we defined an equivalent likelihood discriminant from information on 
the location and penetration depth of the associated track in the KLM. In addition to the 
requirements L e > 0.8 and > 0.95, the tagging lepton tracks were required to satisfy 
1.4 GeV/c < p <2.2 GeV/c, 45° < 6 Uh < 125°, \dr IP \ < 0.2 cm, and \dz IP \ < 10 cm, where 
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p is the CM momentum, 0i a b is the polar angle of the track in the laboratory frame, and 
dzip and drjp are the distances of closest approach to the IP in the direction of the beam 
and in the plane perpendicular to it, respectively. For tagging electrons that satisfy the 
above selection criteria, the detection efficiency, determined by embedding single simulated 
tracks in multi-hadron data, is 92%, and the misidentification probability for hadrons, de- 
termined from Kg — ^ 7t + 7t~ decays, is less than 0.4%. For tagging muons, the efficiency and 
misidentification probability, determined in a similar way, are 86% and 1.2%, respectively. 

To select electron candidates for the spectrum measurement, we required that tracks 
satisfy p^ > 0.5 GeV/c, 46° < 9^ < 125° , \dr IP \ < 0.3 cm and \dz IP \ < 4.0 cm (pi ab is the 
momentum in the laboratory frame). We used dE/dx, TOF, and ACC information to reject 
tracks that are strongly identified as being other than electrons by a likelihood criterion of 
Lnon-e < 0.95. The efficiency of this cut for retaining electrons was 98.6%. We did not use 
ECL cluster information because that was used at a later stage to extract electron yields by 
fitting the E/p distribution. 

Further requirements were made to reduce backgrounds from several known sources. 
Events were rejected if any track satisfying the lepton identification requirement of L e > 0.8 
or L M > 0.95 could be paired with any other oppositely charged track to obtain an invariant 
mass within 50 MeV/c 2 of the J/ip mass. To remove electrons originating from 7 conversions, 
we rejected all pairs of oppositely-charged tracks where at least one track satisfied L e > 0.8 
and the e + e~ invariant mass was less than 0.10 GeV/c 2 . To remove electrons from 7r° Dalitz 
decays, oppositely charged track pairs with e + e~ invariant mass less than 0.10 GeV/c 2 
were combined with photon candidates and rejected if the e + e~7 invariant mass of any 
combination was consistent with the tt° mass. 

To reduce dileptons in continuum events, which tend to be back-to-back, and "ghost" 
tracks (one track mis-reconstructed as two nearly collinear tracks), we required —0.8 < 
cos^e < 0.998 for both opposite-sign and same-sign dileptons, where 0£ e is the opening 
angle between the tagging lepton and the spectrum electron in the CM frame. 

To reject same--B backgrounds, where both leptons originate from just one B meson, we 
look at angular, kinematic, and charge correlations. The spectrum electron and the tagging 
lepton, if correctly selected, have opposite charges unless one of the parent B mesons has 
mixed. If a secondary electron is selected as the spectrum electron, the charge combination 
is opposite to the above case. If a secondary electron from the same B as the tagging lepton 
is selected as the spectrum electron, it has a charge opposite to that of the tagging lepton. 
Opposite-sign dileptons were required to satisfy either p e + cos6>£ e > 1.2 (p e in GeV/c) or 
cos#£ e > 0.3. This criterion reduced the fraction of the s&me-B background in the spectrum 
electron candidates to 2% while retaining 61% of the signal. 

IV. SIGNAL YIELD AND BACKGROUND SUBTRACTION 

To extract the number of dilepton events from the lepton-tagged events, a maximum 
likelihood fit was performed on the E/p distributions of the spectrum electron candidates 
in each momentum (p) and polar angle (#i a b) bin as shown in Fig. [T], where the momentum 
range is divided into 46 bins of width 0.05 GeV/c (0.5-2.8 GeV/c), the polar angle region 
is divided into 4 bins (46-60°, 60-78°, 78-97° and 97-125°), and same- and opposite-sign 
dileptons are treated separately. The signal probability density function (PDF) was obtained 
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using e + and e~ from photon conversions and from the two-photon process. The background 
PDF was taken from the sample of tracks rejected by the electron identification requirements 
for the spectrum study. By determining both signal and background PDF's in the same bins 
of momentum and angle as the dilepton sample, we removed any systematic effects due to 
differences in momentum and angular distribution between the samples used to obtain the 
PDF's and our dilepton sample. 

The raw momentum spectra in the CM frame for both opposite-sign and same-sign 
electrons with the above selections applied are shown in Fig. [| The total yields were 
19722 ± 147 opposite-sign (OS) and 11224 ± 117 same-sign (SS) dileptons in the data taken 
at the T(AS), and 791 ± 105 OS and 448 ± 87 SS below resonance (after scaling). The 
raw yields are shown in Table |. The continuum contribution was subtracted by scaling the 
dilepton yields in the off-resonance data. The scaling factor is the ratio of the on- and off- 
resonance integrated luminosity, after correcting for the energy dependence of the continuum 
cross section. 

These uncorrected spectra contain several backgrounds that are accounted for and sub- 
tracted. Leptons from J/ip decay and electrons from 7 conversion and tt° Dalitz decays 
can survive when one of the pair has escaped detection. These backgrounds are deter- 
mined by Monte Carlo simulation, and the uncertainties were evaluated from the error on 
each rate. Contributions from previously unaccounted sources of secondary leptons from 
B decay were modeled via Monte Carlo, assuming current values of branching fractions: 



B -> Xr+z/,r+ -► e + Y (B(B -► Xt + v) = (2.6 ± 0.4)% @), B -► D S X,D S -> e~Y 
(B(B -> D S X)B(D S -> e~Y) = (0.8 ± 0.3)%), B -> A C X, A c -> e+Y (B(B -> A C X)B(A C -> 
e + Y) = (0.29 ± 0.12)%), B -> DX S ,D -> Ye + v via b -> ccs {B{B -> DX 5 ) = (7.9 ± 
2.2)% False tag leptons arise mainly from secondary leptons and fake muons, whose 

contamination was estimated to be (2.2 ± 0.4)% and (1.5 ± 0.8)%, respectively. All back- 
grounds in the dilepton events were scaled by normalizing to the number of lepton-tagged 
events. The background yields and their errors are summarized in Table |. The spectra of 
the total BB backgrounds are included in Fig. [| 



V. MEASUREMENT OF BRANCHING FRACTION 



The opposite- and same-sign electron spectra after background corrections can be written 



as 



dN+ 



dp 
dN±± 
dp 



efc2(p)^(p)X tag 



dB(B -»■ Xe+v) 
dp 



x dB(B — > DX,D — > Ye~v) 

- "» + — t v -x. 

dB(B^Xe + v) dB(B — > DX, D — >■ Ye~u) 
dp dp 



(1) 
(2) 



where N tag is the number of tag leptons, rj(p) is the electron identification efficiency as 
a function of the momentum, which also includes the momentum smearing effect from 
bremsstrahlung, and €k(p) is the efficiency of our kinematic selection, determined from 
a Monte Carlo simulation. \ is an effective mixing parameter for all B mesons com- 
ing from T(4S) decay, x = Xo/00, where xo is the fraction of neutral B events decay- 
ing as mixed (B°B° or B°B°) and /oo (/+-) is the branching fraction for T(4S) decay 
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into B°B°(B + B-). We used x = 0.0853 ± 0.0055, calculated using X o = 0.174 ± 0.009, 
f + _/f 00 = 1.04_± 0.08 ]T§, and / + _ + / Q0 = 1.0. The two spectra, dB(B -> Xe + v)/dp and 
dB(B — ► DX,D — ► Ye~u)/dp, can be obtained from the measured dN + _/dp and dN±±/dp 
spectra by solving equations (|T|) and (Q) simultaneously. 

The number of tags was determined by counting leptons above 1.4 GeV/c, subtracting 
backgrounds, and correcting for the relative efficiency for selecting dilepton events compared 
to single-lepton events [(94.8±0.5) %]. As a consequence, we found the number of tags to be 
-^Vtag = 528975 ± 1412 (stat). The difference in the relative efficiency between opposite-sign 
and same-sign dilepton events was also corrected for in each electron spectrum. 

By solving the equations, we obtained the momentum spectra for both the primary and 
secondary electrons, as shown in Fig. |3|. The primary electron spectrum was integrated from 
0.6 to 2.8 GeV/c to extract the partial branching fraction for semileptonic B decay, 

B(B -> Xe + u | p > 0.6 GeV/c) = (10.24 ± O.ll(stat) ± 0.46(syst))%. (3) 

The relative size of the undetected region of the electron spectrum was studied using the 
ACCMM [17] and ISGW2 models and a semi-empirical shape [|19j for b — > c decay, 



with options to include semileptonic decay with baryons fl9|| , and charmless semileptonic 



B decays. As part of our evaluation, we fitted the measured primary spectrum to the 
prediction of each model. Results of the fitting are shown in Table [II]. For all models, the 
fits were improved when the proportion of decays B — > D**£u was allowed to float to higher 
values than were nominally set in the model. We obtained the best x 2 fit using a semi- 
empirical shape, with B{B — > X u tv) fixed to 0.167%, and the following quantities allowed 
to vary (fitted values with statistical error are shown): B(B -> LK*W) = 7.31 ± 0.28%, 
B(B -> D**eu) = 2.51 ±1.00%, B(B -> D^heu) = 0.78 ±0.96% and B(B -> baryons ev) = 
0.03 ± 0.31%. We extrapolated the branching fraction below 0.6 GeV/c using the electron 
momentum spectrum of this "best-fit" model. We found that the extrapolated area was 
(6.1 ± 0.7)% of the total decay width, where the error was taken from the largest difference 
between the best-fit model and the other models (excepting ACCMM, which gave a poor 
fit). We obtained: 

B{B -> Xe + u) = (10.90 ± 0.12(stat) ± 0.49(syst))%. (4) 

This result is consistent with other measurements . 

To determine the CKM matrix element \V c b\ we first subtracted the contribution from 
charmless decay (B(B -> XJz/) = (0.167±0.055)%) and assumed B(B -> Xe + u) = B(B -> 
Xlv). Using a formula based on the heavy quark expansion [20|, with jj^ = (0.5±0.2) GeV 2 , 



rrib = 4.58 ± 0.05 GeV/c 2 and the world-average value for B lifetime, tb = 1.607±0.021 ps 
(average of B° and B + ) |l6l , our evaluation is 

\V cb \ = 0.0408 ± 0.0010(exp) ± 0.0025(th), (5) 

where the first error includes both the statistical and systematic uncertainties and the second 
error arises from theoretical uncertainties. 

The branching fraction for secondary decay, BiB —>■ DX, D — > Ye~v), was extracted by 
fitting the secondary electron spectrum to the ACCMM and ISGW2 model predictions. We 
found 
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B(B -> DX, D -> Fe"z>) = (7.06 ± 0.14(stat) ± 0.48(syst))% (ACCMM) (6) 

= (7.84 ± 0.13(stat) ± 0.48(syst))% (ISGW2) 

The obtained branching fraction of the secondary charm semileptonic decay is consistent 
with results from other experiments [|I7^ , |l6ll . 

VI. SYSTEMATIC ERRORS 

The systematic uncertainties on B(B — > Xe + v) and BiB — > DX,D — ► Ye~v) are sum- 



marized in Table |TJ. The main contributions are the uncertainties in the tracking and 
identification efficiencies for the electrons used in the spectrum study. The tracking effi- 
ciency was estimated by embedding electron tracks into hadronic event data and repeating 
the tracking procedures to evaluate the efficiency of the embedded electron. The embed- 
ded electrons were obtained from single electron Monte Carlo as well as from electrons in 
two-photon data. The tracking efficiency was also estimated by embedding into simulated 
BB events, and the absolute difference in values from the two methods was taken as the 
uncertainty. 

The electron identification included two steps: hadron-muon rejection and E/p fitting. 
The inefficiency in the first step (1.4%) was found from a study using photon conversion 
electrons embedded into hadronic events in data. The systematic uncertainty of 0.8% on this 
effect was estimated by taking the difference between results obtained from embedding in 
data and MC. Uncertainties in E/p fitting can come from impurities in the data samples used 
to obtain the signal and background shapes as well as differences in the event environment 
between those data samples and dileptons. The signal shapes were obtained from photon 
conversion electrons in the hadronic sample and two photon data, while the background 
shapes were obtained from tracks failing the hadron-muon rejection. Impurities in these 
samples were estimated via Monte Carlo, and systematic uncertainties of 0.9% and 0.8% 
were assigned for fitting of electron and hadron shapes, respectively. The uncertainty due to 
differences in event environment was also estimated using MC and found to be 1.3%. The 
quadratic sum of the contributions gives a net uncertainty on electron identification of 2.0%. 

The uncertainty in the kinematic selection efficiency was estimated by changing the 
event selection criteria for signal events. Instead of requiring either p e + cos0&. > 1.2 or 
cosOte > 0.3, we required either p e + cos#&. > 1.0 or cos0&. > 0. We also varied the polar 
angle region for tag leptons and considered the ranges; 45° < 9^ < 65°, 65° < 6> lab < 90° 
and 90° < 9^ < 125°. The maximum deviation of the results was assigned to a systematic 
error. 

The systematic uncertainty in the continuum subtraction was attributed to the uncer- 
tainty on the normalization between on-resonance and off-resonance. The uncertainty was 
estimated from the error on the integrated luminosity. This effect is larger for the lepton tag 
events than the dilepton events since the fraction of continuum background for tag events 
is much higher than for dilepton events. 

There were two contributions to the uncertainty on the background subtraction in the 
dilepton events. The first is the uncertainties on the branching fractions of the various 
sources, which are given in Section [TV[ The second is the uncertainty on the scale factor 
used to normalize background yields obtained from Monte Carlo simulation. The dominant 
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contributor to this effect is the uncertainty in the number of tag leptons, which was used to 
determine scaling factors. 

The main sources for the uncertainty on the relative efficiency are the invariant mass 
cuts to reject electrons from J/ip and gamma conversions. This uncertainty was estimated 
by taking the difference of the final yields with and without the cuts. 



VII. CONCLUSION 

In conclusion, we have measured the branching fraction for inclusive semileptonic B 
decay using a high-momentum lepton tag. The result is 

B(B -> Xe + v) = (10.90 ± 0.12 ± 0.49)%, (7) 

where the first error is statistical and the second is systematic. This result agrees with 
other measurements. From this result and the world average B meson lifetime [0 , we have 
extracted the CKM parameter |V^,| as follows: 

\V cb \ = 0.0408 ± 0.0010(exp) ± 0.0025(th). (8) 
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TABLE I. Yield of electrons with momentum of 0.6 GeV/c < p e <2.8 GeV/c in lepton-tagged 
events. The background yields and their statistical errors were determined by a Monte Carlo 
simulation. The second errors are due to branching fraction uncertainties. 





Opposite sign 


Same sign 


On-resonance data 


19722 ± 147 


11224 ± 117 


Scaled off-resonance 


791 ± 105 


448 ± 87 


Continuum subtracted 


18931 ± 181 


10776 ± 146 


Lepton from JM or ib } 

1 ITT 


167 ± 8 ± 


33 


179 ± 8 ± 36 


e from 7 


127 ± 7 ± 


32 


306 ± 10 ± 77 


e from 7r°, 77 


30 ± 3 ± 


8 


88 ± 6 ± 22 


e from r 


406 ± 12 ± 


84 


75 ± 5 ± 16 


e from A c 


8 ± 2 ± 


3 


177 ± 8 ± 71 


e from D s 


293 ± 10 ± 


88 


64 ± 5 ± 19 


e from D(B -> DX) 


285 ± 10 ± 


86 


54 ± 4 ± 16 


Same-f? 


397 ± 12 ± 


79 




Tag from D 


160 ± 8 ± 


32 


600 ± 15 ± 120 


Fake lepton tags 


240 ± 9 ± 120 


323 ±11 ± 162 


Other backgrounds 


43 ± 4 ± 


13 


51 ± 4 ± 15 


Total background 


2156 ± 28 ±215 


1917 ± 26 ± 233 


Background subtracted 


16775 ±183 ±215 


8859 ± 148 ± 233 ± 148± 233 
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TABLE II. Goodness-of-fit for various models fitted to the primary spectrum. Models modified 
to vary the fraction of D** are indicated by '**', and the "best fit" is described in the text. 



Model X 2 /DOF 

QQ98 j@j 50.0/43 

ISGW2 || 70.0/43 

ISGW 48.3/43 

ACCMM 111 125.8/43 



ISGW2** 32.5/41 
QQ98 ** 30.2/42 
QQ98 (best fit) 28.4/39 



TABLE III. Summary of systematic errors on the branching fractions. 



Source of uncertainty 


AB(B ^ Xe+u) (%) 


^-{B -» DX,D -+Ye~v) (%) 


Tracking 


2.9 


2.9 


Electron identification 


2.0 


2.0 


Kinematic selection 


1.4 


4.1 


Continuum subtraction 


0.7 


0.9 


Background subtraction 


1.8 


2.3 


Mixing parameter 


0.6 


1.5 


Relative efficiency 


1.3 


0.8 


Model prediction 


0.7 


10.4 


Total 


4.5 


12.1 
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FIG. 1. E/p distributions for electron candidates in the momentum ranges 1.15-1.2 GeV/c 
(top) and 0.65-0.7 GeV/c (bottom) with the electron candidate track in the polar angle range 
between 97° and 125° in the laboratury frame. The points with error bars are the data, the 
solid histogram is the sum of signal and background components, and the shaded area shows the 
non-electron background contribution. 
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FIG. 2. Electron momentum spectra. The top and bottom figures show the spectra of oppo- 
site-sign and same-sign electrons, respectively. The closed circles are the on-resonance data. The 
triangles show the scaled off-resonance data. The error bars indicate only the statistical error. The 
histogram is the MC-determined BB background. 



14 




p e (GeV/c) 

FIG. 3. Final results for the electron spectra of primary (circles) and secondary (triangles) 
semileptonic decays. The solid curve superimposed on the primary electron spectrum is the best-fit 
model described in the text, while the curve for the secondary electron spectrum is the ISGW2 
model. 
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